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Abstract 

The production of an intermediate-mass Higgs boson in processes of 
the kind e + e~ — > 6 fermions at the energies of future linear colliders 
is studied. The recently developed and fully automatic algorithm/code 
ALPHA is used to compute the tree-level scattering amplitudes for the 
reactions e + e" — ► fi + fi^r^v i T ud, n~e~i> e ud. The code has been in- 
terfaced with the Monte Carlo program HIGGSPV/WWGENPV, properly 
adapted to 6-fermion production, in order to provide realistic results, 
both in the form of cross sections and event samples at the partonic level. 
Phenomenological results, that incorporate the effects of initial-state ra- 
diation and beamstrahlung, are shown and commented, emphasizing the 
potentials of full six-fermion calculations for precise background evalu- 
ation as well as for detailed studies of the fundamental properties of the 
Higgs particle. 
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1 Introduction 



The investigation of the mechanism of electroweak symmetry breaking through 
the search for the Higgs particle constitutes one the main tasks of the experi- 
ments at present and future colliders. Present-day efforts at LEP2 are expected 
to reach the mass limit m h « ^Js — 100 GeV and therefore will be unable to 
pursue the search for the Higgs boson in the interesting range of intermediate 
mass, i.e. Mz < < 2Mz- This mass interval is actually particularly attrac- 
tive since a number of theoretical arguments as well as fits to precision elec- 
troweak data 0, provide indication for a weakly-coupled Higgs boson lighter 
than, say, 200 GeV. On the other hand, the discovery of the Higgs boson in 
the intermediate-mass interval might pose severe problems at the future large 
hadron collider LHC, in spite of continuous progress in detector performances 
and research strategies |§ as well as in updating theoretical predictions [|J. 

From this point of view, the next-generation of e + e~ linear colliders (NLC) 
could be of great help in the complementary search for a Higgs boson with 
intermediate mass, as well as in the determination of its fundamental properties 
and couplings to SM particles. At the NLC, for "moderate" centre of mass 
(cm.) energies up to 0.5 TeV, the main production mechanism for the Higgs 
boson is the Higgs-strahlung process e + e~ — > ZH . When considering the 
M H > 135 -T- 140 GeV range, the above Higgs-strahlung process gives rise to 
a six-fermion final state, originating from Z decaying into // pairs and Higgs 
boson decaying predominantly into WW pairs, with the subsequent decay 
of each W into leptonic lv\ or hadronic ud final states. Therefore, the full 
calculation of the rate of intermediate-mass Higgs boson production at the NLC 
faces with the problem of a full calculation of 2 — > 6 scattering amplitudes, 
including Higgs signal and background processes in the electroweak theory. 
The same difficulty applies to other very important research streams of the 
NLC, namely tt and three-vector bosons production, that both generate six- 
fermion final states. 

A few calculations of such processes have been very recently performed, 
and a short account of the existing approaches can be found in ||. Although 
different in several computational and numerical aspects, all these strategies 
rely upon the calculation of the very large number of Feynman diagrams 
involved (of the order of hundreds), exploiting in particular the experience 
accumulated during the last few years in the computations of the processes 
e + e~ — > 4 fermions developed for W — pair and "light" Higgs boson production 
at LEP2 |], |7|, H H [1(J . Concerning these recent calculations of specific pro- 
cesses of the kind e + e~ — > 6 fermions (6/), in [O the full set of diagrams for the 
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semi-leptonic processes e + e~ — *> e~z/ e 4g and e + e _ — > /i~z/ M 4g have been com- 
puted and physical distributions of interest for top-quark and WW Z physics 
have been analyzed in detail. Similar calculations are in progress following a 
fully computerized approach for the automatic computation of the Feynman 
diagrams by means of the GRACE system |12|]. Previous calculations of helic- 
ity amplitudes relative to multi-particle production processes and relevant for 
collider phenomenology can also be found in the literature ffi 



Here an alternative attack strategy to the problem is presented, using a 
theoretical algorithm (and the corresponding code) recently proposed in the 



literature and known as ALPHA fL4 |. For a given Lagrangian, this algorithm, 
of iterative nature, allows to compute numerically, in a fully automatic way, 
the tree-level scattering amplitudes, without using the Feynman diagrams, 
and it turns out to be particularly powerful for the calculation of those pro- 
cesses involving a high number of final-state particles. The resulting code 
has already been applied with success to the calculation of the rates of multi- 
particle production reactions of interest for the LEP2 and NLC experimental 



programme, such as e e — > 4 fermions, e e — > 4 fermions+7 fll4| , |15| and 



77 — > v e e~ud [16]. In particular, in the case of four-fermion production pro- 
cesses the numerical results obtained with ALPHA have been compared in de- 
tail with those of independent formulations based on standard computational 
techniques, showing excellent agreement and thus providing a stringent test 
of the algorithm and of the resulting code [J], [J. In order to produce realis- 
tic results, the algorithm has been interfaced with the Monte Carlo program 
HIGGSPV/WWGENPV [TP?], |T8[, developed for WW and Higgs boson physics in pro- 



cesses with four fermions in the final state in the context of LEP2 physics [[19] 



and properly adapted to 6/ production. The computational tool developed 
can work both as an integrator of weighted events, providing cross sections 
for any given experimental set up, and as a generator of unweighted events, 
providing event samples suited for physics analysis and/or detector simulation. 

To the best of the authors' knowledge, a full six-fermion calculation of 
intermediate-mass higgs boson production in e + e~ collisions and relative phe- 
nomenological analysis have not yet been performed. In this paper the results 
obtained for two different channels are presented, trying to emphasize the 
usefulness of a complete 6f calculation from the point of view of precise back- 
ground evaluations in Higgs boson searches as well as for the determination of 
the fundamental properties (mass, spin, etc.) of this particle. Recent reviews 
on the main aspects of Higgs boson phenomenology at future e + e~ linear col- 



liders can be found in [2CJ, |2J, |22|, |23|, [24], |2q| . Moreover, also full calculations 
for processes of the kind e + e~ — *> bbWW, bbZZ and Ajets + W, of interest for 
Higgs boson searches at NLC, but with on shell final-state vector bosons, can 
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be found in [E6l E7 



The paper is organized as follows. In Section 2, the details of the calcu- 
lation are presented, considering intermediate-mass Higgs boson production 
in the six-fermion reactions e + e~ — > fi + fi~T~v T ud,n + [i~e~v e ud; a sample of 
illustrative numerical results is shown and commented in Section 3 and the 
main conclusions as well as possible perspectives are drawn in Section 4. 



2 Six-fermion Calculation 

In order to study the intermediate-mass Higgs boson production at the NLC 
in channels that are as free as possible from large background contamination, 
only the processes 

e + e~ — > ji + fi~T~v T ud, /i + \T e~v e ud (1) 

are considered, which differ for their content of electroweak backgrounds but 
are both completely unaffected by QCD backgrounds. 

In the following, some details concerning the strategies followed for the 
calculation of the physical amplitudes and the phase space of the above 2^6/ 
processes will be described. 

As said in the introduction, to compute the matrix elements, a recently 
proposed technique (ALPHA) has been employed. It doesn't require the eval- 
uation of Feynman graphs and allows to obtain automatically, according to 
an iterative procedure, the tree-level scattering amplitudes of any given la- 
grangian. This is achieved by exploiting the relation between the generator 
of one-particle irreducible Green functions and the generator of the connected 



Green functions. The interested reader is referred to the original literature [14 
for a detailed description of the method. Here it is worth pointing out that 

• the procedure, which is entirely automatic, has been checked successfully 
against the precise predictions existing for the processes e + e~ — > 4/ |7|, |9| 
and applied to obtain original results for the rate of the reactions e + e _ — > 
4/ + 7 [ 15 1 and 77 — > 4 fermions [HJ; 



to calculate the amplitude for a certain process, an input file is provided, 
where one has simply to specify the type of the process and the total 
number of particles involved in intial and final states. For the pure QED 
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case the required input is an integer vector N prt = (m, . . . ,n 6 ) where 
ni,...,n,Q are the numbers of initial e~, final e + , final e~ , initial e + , 
initial and final photons respectively; once this vector is initialized and 
the proper kinematics is provided ALPHA returns the scattering matrix 
element s;[] 

• the present calculation is the first application of ALPHA to the compu- 
tation of e + e~ — > 6f processes; this enables, in particular, to check yet 
untested parts of the input lagrangian (i.e. the coupling of the higgs bo- 
son to the W and Z bosons and some of the quartic non-abelian gauge 
boson couplings) and to show the potentials of the algorithm for multi- 
particle production processes involving a very large number of diagrams. 

The kinematics of the 2 — > 6/ processes has been treated generating the 6- 
body phase space recursively. This choice is particularly useful since processes 
where a particle decays into other particles which subsequently decay are being 
considered (see fig. [I]). In particular, denoting by the four-momentum 

of the incoming electron (positron), q i} i = 1, ... ,6, the final-state outgoing 
four- momenta and P = p_ + p + , the phase space can be written as 

rf$ 6 (P; qi , . . . , q 6 ) = (2vr) 12 d$ 2 (P; Q z , Q h )d$ 2 {Q z ; q u q 2 ) 

d$ 2 (Qh; Qw+,Qw-)d$2{Qw+; <?3, 94)^2 (Qw-; <?5, <?e) 
dQ 2 z dQ 2 h dQ 2 w+ dQ 2 w ., (2) 

where Q z = <?i + <?2, Qh = Ei=3 Qh Qw+ =Q3 + q4 and Q w - = q 5 + q 6 . Globally, 
fourteen independent variables (including a trivial overall azimuthal angle) are 
required; they have been chosen according to the following scheme: 

• four invariant masses Q z , Q\, Qw+j Q\y- 

• five and <^-angle pairs, in the rest frame of each decaying "particle" , 
namely in the frames where P — 0, Q z = 0, Qh = 0, Qw+ — and 
Qw- — 0, respectively; once these variables are generated, the four- 
momenta of the outgoing fermions are derived in each "rest" frame and 
then boosted back to the laboratory frame. 

The choice of the kinematics has been motivated by the dynamics of the Higgs 
boson signal; the same phase space is also conveniently employed for the back- 
ground processes. 

1 For the standard model the only difference is that the vector N prt has much more entries 
(a pair for each SM particle). 
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In order to obtain reliable phenomenological results for the energies of the 
future e + e~ colliders, the lowest-order calculation of the matrix elements has 
to be supplemented with the potentially large effects of energy losses due to 
the emission of initial state photons, caused by the mutual interaction of the 
colliding electrons and positrons at the "parton" level (bremsstrahlung or ISR), 
as well as to the synchrotron radiation generated by the strong electromagnetic 
interaction between dense bunches of particles (beamstrahlung). Both the 
contribution of ISR and beamstrahlung have been incorporated in the physics 
simulation in order to provide as realistic as possible predictions. The task 
has been accomplished by interfacing ALPHA with the Monte Carlo program 
HIGGSPV/WWGENPV, properly adapted to treat 6/ production processes. This 
allowed to exploit as much as possible the expertise developed in the context 
of the LEP2 workshop ||19|| , concerning 4-fermion production processes. In 



particular, the method of QED structure functions (SF) [ES| is used to account 
for ISR according to the factorized formula 

a ( s ) = E / Q dxi dx2 D ( Xl ' S ) D ( X *> s ) d [ PS ] rfjpsjyp' ( 3 ) 

where D(x,s) is the electron SF, d[PS] denotes the volume element of the 
6-body phase space, da/d[PS] is the tree-level hard scattering density in the 
cm. frame, built by means of the total scattering amplitude as returned by 
ALPHA and of the proper jacobian factors; Wi and W — J2i w i are weights in- 
troduced to increase the efficiency of the Monte Carlo integration, namely by 
incorporating the weights Wi in the integration measure according to a multi- 
channel importance sampling PDfl; the random number generator employed is 



RANLUX [30]. Beamstrahlung corrections are simulated using the parameteri- 



zations of such effects recently implemented in the library CIRCE pl| , namely 
according to the formula 



<*) = E 

i 



1 rl 

dzidz 2 D BS (zi 1 z 2 ; y/s) / dx\dx 2 D(x\, s)D(x 2 , s)d[PS]- 

Jo i 



da Wa 



d[PS]W 
(4) 

where Dbs is the beamstrahlung distribution.^ When ISR and/or beam- 
strahlung are taken into account, the four-momenta of the final-state particles 
are additionally boosted to the laboratory frame, in order to take into account 
the cm. boost generated by such effects. At present, ISR is treated in the 
strictly collinear approximation. 

2 The TESLA accelerator code has been used. The library CIRCE provides the beam- 
strahlung parameterization for some fixed y/s, in particular y/s = 350 and 500 GeV. For 
one of the cm. energies considered in the following simulations, namely y/s = 360 GeV, 
the parameterization relative to 350 GeV has been used, as explicitly allowed in the original 
documentation. 
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A technical but important item in the numerical calculation is the choice 
of the weight functions Wi employed for the implementation of the importance 
sampling. This choice has to be driven by the physics involved. When con- 
sidering only the Higgs boson signal (see fig. 0), it has to be observed that 
muon pairs come from the Z boson decay, while the remaining four final-state 
fermions come from the Higgs boson decay. Moreover, the hadronic pair comes 
from the W + boson decay and the additional leptonic pair from the W~ bo- 
son decay. This means that the most populated phase space region is in the 
surroundings of Q\ ~ Mf , Q\ ~ m\ and Q\ v ± — M^. When considering 
the leading background processes common to both the channels examined (see 
for instance fig. ^j, the muon pairs can come also, for instance, from photon 
conversion; moreover, Q\ is no more peaked around m\\ the hadronic and the 
additional leptonic pairs still come from Ws decay. At last, when considering 
background processes typical of the channel with the electron in the final state 
(see for instance fig. |j) the ev leptonic pair invariant mass is no more peaked 
around M^. Following these considerations, the importance sampling adopted 
has the following features: 

• Q\ is sampled piecewise, according to a Breit-Wigner (BW) density 
around Q\ = ml and a flat density elsewhere; the relative weights of 
the BW and flat densities are tuned according to the signal/background 
ratio for any given y/s; 

• Q 2 Z is sampled piecewise, according to a BW density around Q 2 Z = Mf , 
according to 1/Q% in the lower tail and a flat density in the upper one; 

• Q\ v ± are always sampled according to a BW density; 

• since a realistic event selection requires a minimum scattering angle for 
the leptons, no particular sampling has been performed in the final- 
electron forward region; in particular, a minimum angle of 5° between 
the final electron and the incoming beams has been imposed, in such a 
way that the t-channel "singularity" (see for instance fig. |]) is excluded. 

The "soft-photon singularity" of the electron SF's is sampled in the stan- 



dard way, as can be found in the literature [17, IS 



CPU performances have not been the main concern of the study. Anyway, 
it is worth saying that the Monte Carlo code developed, in the worst case, 
generates around 3500 unweighted events in 8h CPU and, as integrator of 
weighted events, provides a cross section with a statistical error of around 1% 
in 3h CPU on a DEC ALPHA station. Such performances are sufficient for 
any realistic simulation. 
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3 Numerical results and discussion 



In order to produce the numerical results shown in the present section, the 
input parameters used are the Fermi coupling constant Gp, the Z boson mass 
M z = 91.1888 GeV and the W boson mass M w = 80.23 GeV. The weak 
mixing angle and the QED coupling constant are then computed as sin 2 $ w = 
1 - M^/Mf and Ana = g 2 sin 2 $ w , where g 2 = 8M^G F /y/2. Given these 
input parameters, all the other relevant quantities, such as the Z and W 
bosons widths, have been computed at the tree-level approximation and the 
massive boson propagators are chosen as ~ 1/ (p 2 — M 2 + iTM ) . This choice of 
input parameters has been performed in order to avoid problems concerning 
violations of the Ward identities, or, more precisely, to confine all possible 
gauge violation effects into the inclusion of the finite fixed widths for the 
massive bosons. 

As far as the Higgs boson width is concerned, in the present simulation the 
following contributions have been taken into account: the fermionic contribu- 
tions h — ► fifi, tt, cc, bb, where the hadronic widths have been QCD corrected 
taking into account the running quark mass effects ||; the gluonic contribution 
h — > gg according to ||; the vector boson contribution h — > V*V*, according 
to ref. 131 



The result has been tested to be U(l) em gauge- invariant, by comparing the 
results obtained with two different inverse photon propagators in the effective 
lagrangian implemented in ALPHA. A priori, one could expect £7(1) gauge in- 
variance problems to be severe in particular in the channel with the electron 
in the final state. Anyway, such potential problems |33|] are avoided requiring 



a realistic event selection, as described above. Concerning the SU{2) gauge 



invariance, it has been pointed out in [34] that the fixed- width scheme, used 



in the present study, is not SU(2) gauge invariant and could in principle be 
responsible of a bad high-energy behaviour in processes involving six fermions 
in the final state. The numerical relevance of the SU(2) gauge violations of 
the (fixed widths) results has then been checked by comparing them with re- 
sults obtained in the so called "fudge scheme" f35|, that is by construction 



gauge-invariant, and finding them compatible at the per cent level. 

All the numerical results, but those shown in fig. [5], have been obtained 
by imposing the following a priori cuts: the invariant mass of the muon pair 
larger than 20 GeV; the invariant mass of the hadronic system larger than 
10 GeV; the angle between the charged leptons and the beams larger than 5°. 

The results shown in figs. [5] and ^| have been obtained by using the program 
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Monte Carlo integrator of weighted events. 



Figure ^| shows the comparison between the fully extrapolated cross sections 
for the Higgs boson signal alone, as obtained by means of the present complete 
6f calculation on the one hand, and within the narrow-width approximation 



(NWA) according to ref. [36| on the other one. The two windows show the 
comparison as a function of the Higgs boson mass, for two fixed ^fs = 360 
and 500 GeV, and as a function of the cm. energy, for three fixed Higgs 
boson masses rrih = 150, 160 and 170 GeV, respectively. In both cases the 
relative deviation R = cr^f / &nwa — 1 is considered. In particular, in the first 
plot two cases are considered, namely the case in which the total Higgs boson 
decay width has been used in the full 6f calculation (solid and dashed lines), 
and the case in which the Higgs boson width in the full 6f calculation is 
identified with T(h — > WW*) (dotted and sparsely dotted lines). The dotted 
and sparsely dotted lines approach zero in the region of "small" Higgs boson 
masses: this represents a check that the 6/ calculation reproduces the results 
already present in the literature for the NWA. The solid and dashed lines 
do not approach zero, because of the presence of the fermionic and gluonic 
contributions to the total Higgs boson width, but this is the realistic case. 
Anyway, the off-shellness effects are contained within a few per cent, as can 
be also seen in the second plot. 

Figure ^ shows the comparison, at the tree-level approximation, between 
the full cross section, the separate contribution of the Higgs boson signal 
and the background only, for two values of the Higgs boson mass = 
150, 170 GeV and as a function of the cm. energy. The first plot shows 
the results for the channel e + e~ — ■> [i + [i~T~v T ud, the second one for the chan- 
nel e + e~ — > fi + fi~e~u e ud. As a first comment, it has to be noticed that, for 
the event selection considered, the channel containing the electron in the final 
state is not very much different from the one containing the r lepton; this 
in spite of the fact that the electron channel receives contributions also from 
t-channel processes of the kind shown in fig. |]. The t-channel background 
processes start to become visible only in the high-energy tail of the plot, say 
around y/s ~ 500 GeV. The statistical error due to the numerical integration 
is of the order of 0.5% and hence invisible for the scale adopted in the plots. 
In both channels, the full calculation is consistent with the incoherent sum of 
signal and background, at the statistical level considered. Signal/background 
interference effects can become relevant for resolutions better than 1%. It has 
to be noticed that, in both cases, yfs ~ 500 GeV represents the turning point, 
after which the background processes become larger than the signal for the cuts 
considered. It has been checked that, in the channel e + e~ — > fi + ^r^u T ud, 
the high-energy behaviour of the cross section is in good agreement with the 
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one found in ref. |37|. For y/s = 360 and 500 GeV also the combined effect 



of ISR and beamstrahlung has been studied. For the full signal+background 
inclusive cross sections, radiative effects turn out to be confined within a few 
per cent, being dominated by ISR. This can be easily understood, since the 
Born approximation inclusive cross section is almost flat, as a function of y/s, 
above, say, 300 GeV. 

The results shown in figs. correspond to y/s = 360 GeV, except for 



figs. H and [13|, where y/s = 500 GeV. They have been obtained by analyzing 
the six-fermion unweighted event samples simulated by using the program as 
a generator. In all the cases, the number of events shown is normalized to 
a fixed integrated luminosity, in order to perform a comparison consistent 
with the corresponding cross sections. The channel e + e~ —>■ /x 4 "fi~ ~ev e ud has 
been considered. As general comments, in these figures the dotted histograms 
represent the Born approximation results, the dashed ones represent the results 
including the effect of ISR and the solid ones represent the full prediction, 
including also the effect of beamstrahlung. It has to be noticed that, among 
the various distributions shown, the only ones that are very sensitive to the 
effect of ISR and/or beamstrahlung are the ones analyzed in figs. [7] and |8]. 

Figure [T| shows the distribution of the missing mass, defined as M miss = 
(P — Qz) 2 , P being the total incoming nominal four-momentum, i.e. before 
the energy losses due to ISR and/or beamstrahlung. In the absence of radiative 
effects, the missing mass coincides with the invariant mass of the Higgs boson 
decay products. This is an experimentally relevant quantity, since the Higgs 
boson width is too tiny to be resolved directly. The radiative effects amount 
generically to populate the high-mass tails of the distribution. A comparison 
between figs. [F] and |8| shows that such effect is larger at y/s = 500 GeV. 

In fig. ^ the invariant mass of the \i pair system is shown, in a window of 
some Z boson widths around the Z boson mass. The distribution is always 
peaked at the Z boson mass; the presence of the Higgs boson signal is revealed 
by the height of the peak, which is enhanced by a factor of about 6 with respect 
to the background. 



Figure [10] shows the distribution of the angle between the total /i + yU~ three- 
momentum and the beam. In the presence of the Higgs boson signal only, it 
coincides with the Z boson scattering angle. As the Higgs boson mass varies, 
the distributions are quite similar, all smoothly peaked around = n/2. The 
presence of the Higgs boson signal is revealed by the height of the peak also 
in this case. 

In fig. [n] the distribution of the variable £i = (cos$* u -|-cos$* d )/2 is shown. 
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The angles are denned as follows. First, the u e "three-momentum" is recon- 
structed as vq = — Yh=i Qi an d its energy component is defined by means of 
the shell relation v$ = |u 6 |. Vq coincides with the true v e four-momentum only 
in the Born approximation; in the presence of radiative effects in the collinear 
approximation Vq is the total lost three-momentum. Then the four vector h 
defined as h = q 3 + g 4 + q 5 + v 6 is considered; in the Born approximation and in 
the presence of the Higgs boson signal only, it coincides with the Higgs boson 
four-momentum. The angles $* are measured in the reference frame in which 
h = 0; in the Born approximation and in the presence of the Higgs boson signal 
alone, it is the rest frame of the Higgs particle. In the presence of the Higgs 
boson signal, the distribution of £1 sharpens around zero. In the above proce- 
dure, there is an obvious arbitrariness in defining actually, it could also be 
defined by means of the energy-momentum conservation, in such a way that 
t>6 is the total lost four momentum. It has been checked that the results are 
not critical with respect to this choice. It has to be noticed that the variable 
£1 is a particularly sensitive variable; in fact, the presence of the Higgs boson 
signal is revealed by a significant change in the shape of the distribution. 



Figures [12] and [13] show the distributions of the variable £2 defined in an 

analogous way to £1 as £ 2 = q%-Qhad/<lT /Qhadi where Q5 and QLd are the four - 
momenta of the outgoing electron and hadronic system, respectively, in the 
reference frame defined above. Incidentally, Qhad coincides with the W + boson 
four-momentum. The cm. energies considered are y/s = 360 and 500 GeV. 
For both the cm. energies, the variable £2 shows an enhanced sensitivity to the 
presence of the Higgs boson particle. Actually, the background distribution 
has a more pronounced peak in the region £2 — — 1- It has to be noticed 
that at y/s = 500 GeV there is a two-peak structure, where the peak around 
£2 — is due to the Higgs boson signal, whereas the one at £ 2 — — 1 is due to 
the background processes, which in this case give a contribution to the cross 
section comparable to the signal contribution (see fig. |6|). 

The shapes of the variables £1,2 are very sensitive to the presence of the 
Higgs boson signal. However, they are only partly determined by the under- 
lying dynamics: in fact they are not sensitive to the spinless nature of the 
Higgs particle. In order to point out the spinless nature of the Higgs boson, 
it is necessary to investigate more exclusive forms of angular correlations. To 
this aim, the variable £ 3 has been defined as cosi?* u or cos$* d , i.e. the co- 
sine of the angle between the charged lepton and the up- or down-quark in 
the reference frame defined above. Figures [TJ] and |I5| show the distribution of 
the two cosines for y/s = 360 GeV and 500 GeV, respectively. In all the his- 
tograms the effect of ISR and beamstrahlung have been included. The shaded 
and hatched histograms represent cos$* u and cos$* d , respectively. The white 
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histogram represents the sum of the previous histograms. In all the cases, 
cos$* u is peaked at —1, both in the presence and in the absence of the Higgs 
boson signal. On the contrary, cos$* d is peaked at —1 in the presence of the 
background processes only, whereas in the presence of the Higgs boson signal 
the region cos $* d — 1 is significantly populated. The shapes of the variables £3 
are sensitive to the spinless nature of the Higgs particle. Such a property has 
been pointed out also in a recent paper as a useful tool in the search for 
the Higgs boson in the decay channel h — > W + W~ — > l^l'^vivy at the LHC. 
Anyway, in that paper the study has been performed on signal and background 
separately, without relying upon a full multi-particle calculation. It is worth 
noting that the sensitivity of the shape is definitely larger at sfs = 360 GeV 
rather than at yfs = 500 GeV, where signal and background are of the same 
order (see fig. H). 



4 Conclusions and outlook 

The production of a Standard Model Higgs boson in e + e~ collisions at very 
high energies has been studied, considering the case of intermediate-mass Higgs 
boson which yields a six-fermion final state. The calculation of the Higgs bo- 
son signal as well as of the background matrix elements has been carried out 
using the automatic algorithm ALPHA, which reveals to be particularly powerful 
for such six-fermion calculations involving in the conventional diagrammatic 
approach a huge number of graphs. The tree-level scattering amplitudes re- 
turned by ALPHA have been corrected for the effects of bremsstrahlung and 
beamstrahlung, in order to provide as comprehensive as possible results, by 
interfacing ALPHA with the Monte Carlo program HIGGSPV/WWGENPV, properly 
adapted to treat 6f production processes. 

As a summary of the main outcomes of the study, the following items have 
to be pointed out: 

• the full 6f calculation as compared with the calculation carried out 
within the narrow- width approximation shows that the various off-shellness 
effects can be of the order of some per cent; for experimental accuracies 
better that some per cent the narrow-width approximation becomes in- 
adequate; 

• the full calculation as compared to the incoherent sum of signal and 
background reveals that, for the channels considered, signal/background 
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interference effects are compatible with zero at the level of precision of 
about 0.5%; 



the missing mass distribution, relevant for the determination of the Higgs 
boson mass, is strongly affected by initial-state radiation and/or beam- 
strahlung effects (radiative effects), as expected; 

the shapes of all the other distributions considered in the paper are only 
smoothly affected by radiative effects; 

the variables that are mostly sensitive to the presence of the Higgs par- 
ticle decaying into a VT-boson pair are £1 and £ 2 of figs. [11], [12] and [13]; 
the variable most sensitive to the spinless nature of the Higgs boson is £3 



of figs. [TJ and [15]; in all the cases, they are variables concerning the spin 
correlations of the Higgs boson decay products, and can be meaningfully 
analyzed only by means of a complete 6/ calculation. 



As a conclusion, the approach described in the present paper represents 
a valuable tool for Higgs boson searches at NLC and, more generally, for the 
study of scattering processes involving multi-particle production. 
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Figure 1: The decay chain in the channel e + e — > // + // -urfr (e )kV(e)- 
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Figure 5: The comparison between the results of the full six fermion simulation 
and the results obtained in the narrow-width approximation, as function of the 
Higgs boson mass m h and of the cm. energy. 
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Figure 6: The cross sections for the channels e + e — > r u T ud and e + e 
fi + fi~ e~ v e ud as functions of the cm. energy. 
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Vs - 360 GeV 




Figure 7: The missing mass distributions at y/s = 360 GeV. Dotted, dashed 
and continuous histograms correspond to Born approximation, Born+ISR and 
Born+ISR+beamstrahlung results, respectively. 
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gure 8: The missing mass distributions at ^Js = 500 GeV. Notation as in 
•0. 
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Figure 9: The /i-pair invariant mass distribution in the surroundings of 
M(/i + //~) = M z . Notation as in fig. |7|. 
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Figure 10: The distribution of the angle between the total /i-pair three- 
momentum and the beams. Notation as in fig. |7|. 
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= 360 GeV 




Figure 11: The distribution of the variable £1 (see the text for the definition) 
at y/s = 360 GeV. Notation as in fig. [7|. 
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Va - 360 GeV 




Figure 12: The distribution of the variable £2 (see the text for the definition) 
at y/s = 360 GeV. Notation as in fig. [7|. 
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Figure 14: The distribution of the variable £ 3 (see the text for the definition) 
at s/s = 360 GeV. All the histograms include the effect of ISR and beam- 
strahlung. The shaded and hatched histograms represent cos-#* u and cos$* d ; 
the white histogram is the sum of the previous ones. 
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